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ABSTRACT 

We present new proper motion and parallax measurements obtained with the European VLBI Network (EVN) at 5 GHz 
for the three isolated pulsars B1929+10, B2020+28, and B2021+51. For B1929+10 we combined our data with earlier 
VLBI measurements and confirm the robustness of the astrometric parameters of this pulsar. For pulsars B2020+28 
and B2021+51 our observations indicate that both stars are almost a factor of two closer to the solar system than 
previously thought, placing them at a distance of 1 .39*Qgg and 1.251 q' 44 kpc. Using our new astrometry, we simulated 
the orbits of all three pulsars in the Galactic potential with the aim to confirm or reject previously proposed birth 
locations. Our observations ultimately rule out a claimed binary origin of B1929+10 and the runaway star £ Ophiuchi 
in Upper Scorpius. A putative common binary origin of B2020+28 and B2021+51 in the Cygnus Superbubble is also 
very unlikely. 

Key words, pulsars: individual: B1929+10, pulsars: individual: B2020+28,pulsars: individual: B2021+51, proper mo¬ 
tions, parallaxes, techniques: interferometric 


1. Introduction 

Typical transverse velocities of i solated pulsars are of the 
order of several hundred kms" 1 (ICordes fc Chernofll ll998: 
Arzoumanian ct al. l2002t iHobbs et al.ll2005lb while those 
of their progenitor O- and B-stars are at most several tens 
of kms^ 1 . In the standard neutron star formation scenario 
this discrepancy is explained by an asymmetry in the super¬ 
nova (SN) explosion that imparts a kick to the forming cen¬ 
tral com pact object, accele rating it to the observed veloci¬ 
ties (e.g. lScheck et al.ll200fii) . As a result of the short lifetime 
of SN remnants (< 10 5 yr) and the typical characteristic 
age of young pulsars (r c ~ 10 6 ”' Myr), direct associations 
between SN-remnants and pulsars are rare. Measurements 
of accurate proper motions p and parallaxes n of pulsars 
can, however, indicate the birth locations of pulsars. The 
combination of both p and the distance d, = 1 /7r yields the 
physical transverse velocity, Vj_, which, given an estimate 
of the radial velocity, V r , allows calculating a trajectory 
that traces the pulsar back to its possible birth location. 
Hence, kinematic ages - as opposed to characteristic ages 
t c = P/2P - of pulsars can be determined and conclusions 
about neutron star formation scenarios can be drawn. 


* franz.kirsten@curtin.edu.au, Member of the International 
Max Planck Research School (IMPRS) for Astronomy and As¬ 
trophysics at the Universities of Bonn and Cologne 


One of the first to calculate pulsar orbits was IWrightl 
(11979D . claiming that the pulsars B1929+10 (J1932+1059) 
and B19 52+29 originate from a fo rmer binary system. More 
recently, iHoogerwerf et al. 1 (|200lf ) used the 3D space veloc¬ 
ity of high-velocity runaway stars and parallax and proper 
motion measurements of young nearby pulsars to extrapo¬ 
late their trajectories back in time. Their simulations indi¬ 
cated that the runaway O-star C, Ophiuchi (£ Oph, HIP 
81377) and the young pulsar B1929+10 were likely to 
have been in a binary system in Upper Scorpius (Scorpius- 
Centaurus association) until about 1 Myr ago. According 
to their analysis, the system was disrupted when the pro¬ 
genitor star underwent a supernova explosion. During that 
event, the space velocity vectors of both ( Oph and the pul¬ 
sar were modified to point away from Upper Scorpius. The 
parameter range for which such a scenario is possible is, 
however, rather small. Improved measurements of p and 7r 
for B1929+10 obtained with the NRAO Very Long Baseline 
Array (VLBA) led to the conclusion that a co mmon origi n 
of the pulsar and C, Oph is unUkehOjChatterjee etalj 2004H . 
Adopting the measurements of lChatteriee et all (12004! ), but 
increasing the reported un certaintie s by factors between 
10 and 30, iBobvlevI (|2008f) and also Tetzlaff et al] (1201011 
repeated the simulations of IHoogerwerf et alT i 200 lh . re- 
postulating a binary origin of B1929+10 and ( Oph in Up¬ 
per Scorpius. 
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Table 1: Observing epochs and arrays. 


Date(s) UT range Array 0 


2010 May 28-29 

2130-0930 

Ef 

Wb 

Jb 

Jv 

On 

Me 

Tr 



Zc 

Bd 

Ur 

Sh 

2010 Oct 27 

1130-2330 

Ef 

Wb 



On 

Me 

Tr 

Ys 

Sv 

Zc 

Bd 



2011 Mar 09 

0300-1500 

Ef 

Wb 

Jb 

Jv 

On 

Me 

Tr 

Ys 

Sv 

Zc 

Bd 

Ur 

Sh 

2011 Jun 01-02 

2120-0920 

Ef 

Wb 

Jb 


On 


Tr 

Ys 

Sv 

Zc 

Bd 

Ur 

Sh 


Notes. O) Ef = Effelsberg, DE, 100 m; Wb = Westerbork Synthesis Radio Telescope, NL, 12 — 14 x 25 m; Jb = Jodrell Bank 
Lovell Telescope, UK, 76 m; Jv = Jodrell Bank Mark2 Telescope, UK, 25 m; On = Onsala, SE, 25 m; Me = Medicina, IT, 32 m; 
Tr = Torun, PL, 32 m; Ys = Yebes, ES, 40 m; Sv = Svetloe, RU, 32 m; Zc = Zelenchukskaya, RU, 32 m; Bd = Badary, RU, 32 m; 
Ur = Urumqi, CN, 32 m; Sh = Shanghai, CN, 25 m. 


Table 2: Calibrator details for each pulsar. 



Pointing centre 

Distance to 

Flux density 6 

Source 

RA 

Dec 

pulsar [deg] 

[mjy beam” 1 ] 

B1929+10 

19:32:14.0160 

10:59:32.868 


50 

J1928+0848 0 

19:28:40.8555 

08:48:48.413 

2.35 

160 

J1934+1043 

19:34:35.0256 

10:43:40.366 

0.63 

50 

B2020+28 

20:22:37.0697 

28:54:22.976 


30 

J2020+2826 0 

20:20:45.8707 

28:26:59.195 

0.61 

70 

J2023+3153 

20:23:19.0173 

31:53:02.306 

2.98 

900 

B2021+51 

20:22:49.8596 

51:54:50.400 


80 

.12025+5028° 

20:25:24.9725 

50:28:39.536 

1.49 

110 

J2023+5427 

20:23:55.8440 

54:27:35.829 

2.55 

500 


Notes. U) These are the calibrators referred to as primary calibrators in the text. For the pulsars this is the apparent pulsar 
flux density measured by employing pulse gating. 


In a similar investigation. IVlemmings et al. (1200 41 iden¬ 
tified the two pulsars B2020+28 and B2021+51 as candi¬ 
dates for a common-origin scenario based on proper mo¬ 
tion and pa rallax measuremen ts obtained with the VLBA 
at 1.4 GHz dBrisken et al.l[2002i) . The authors simulated the 
trajectories of the two pulsars back in time and concluded 
that they most likely originated from a binary system in the 
Cygnus Superbubble that was disrupted when the younger 
of the tw o pulsars was b orn in a supernova. The measure¬ 
ments bv iBrisken et al.1 (1200211 are, however, based on five 
observations covering a time span of only roughly one year. 

Here, we present new measurements of p and n for 
the three pulsars B1929+10, B2020+28, and B2021+51 ob¬ 
tained with the European VLBI Network (EVN) at an ob¬ 
serving frequency of 5 GHz. These observations extend the 
time baseline to more than ten years, allowing for an ex¬ 
tremely high precision in measurements of p and ir for all 
three pulsars. We used these data and ran new simulations 
of trajectories to shed new light on the proposed binary 
origin of B1929+10/C Oph and B2020+28/B2021+51. 

2. Observations and data reduction 

The observations described here were conducted with the 
EVN under project code EV018(A-D). All observations 
used a frequency range 4926.49 — 5054.49 MHz with dual 
circular polarizations and two-bit sampling, for a total 
bit-rate of 1 Gbps per station. The 128 MHz frequency 
range in each polarization was split into eight 16 MHz 
baseband channels. We conducted four epochs of observa¬ 
tions between May 2010 and June 2011, as summarized 
in Table [T] This table also lists the EVN stations that 


successfully participated in the array at each epoch. 

We observed each of the three pulsars in all four 12- 
hour epochs, using phase-referencing. Table [2] summarizes 
the pulsars and characteristics of their phase-reference cali¬ 
brator sources. Our phase-referencing tactics included (i) a 
basic four-minute cycle alternating between the target pul¬ 
sar (2.5 min) and the primary calibrator (1.5 min), and (ii) 
insertion of an additional 1.5 minute scan of the secondary 
calibrator in every second cycle. For bandpass calibration, 
we observed the quasars J1800+3848 and 3C454.3: the 
former about two hours from the beginning and the latter 
four hours from the end of each epoch. This observing 
pattern provides about two hours of integration on each 
pulsar per epoch, yielding a nominally expected sensi¬ 
tivity of the arrays ranging from 10.9 to 14.8 p Jy per beam. 

The data from the telescopes were correlated on the 
EVN software cor relator at the Joint In stitute for VLBI 
in Europe (SFXC; iKeimpema et al.ll2015l) . Each of the 16 
MHz baseband channels was correlated with 32 frequency 
points and one-second coherent integrations. Pulsar scans 
were correlated using the gating/binning capability of 
SFXC. Because these were the first observations to employ 
this mode, we devote a few sentences to describe it. 
Given an ephemeris of a pulsar, SFXC can apply a gate, 
defined by a start and stop fraction of a period, such that 
correlation accumulates only during the in-gate interval. 
Before gating, the pulsar data are de-dispersed. Here, we 
used incoherent de-dispersion (a constant correction per 
0.5 ~MHz frequency point). Coherent de-dispersion has 
subsequently been developed on SFXC. The gate itself 
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Fig. 1: Astrometric accuracy as a function of angular sepa¬ 
ration between phase cali brator source and target . This is a 
reproduction of Fig. 3 in lChatteriee et ahl (j2004f ) (referred 
to as C04 here) to which we added data from our obser¬ 
vations (open and filled stars). Open symbols denote the 
median scatter about the average position of the primary 
phase calibrator obtained from calibrating its visibilities 
with solutions from the secondary calibrator. Filled symbols 
show the median scatter of the observed positions about the 
best-fit model for proper moti on and parallax. The 1.6 GHz 
data (squares) ar e taken from IChatteriee et al.l (|2001h and 
IVlemmings et all (120031) . while the previous 5 GHz data 
(circles) are taken from lGhatteriee et al.l (1200411 . The solid 
line is a weighted least-squares power-law fit to all 5 GHz 
data points. The almost linear relationship between astro- 
metric error and calibrator throw seems to break down be¬ 
yond an angular separation of more than about four de¬ 
grees. 

may be divided further into a number of equal-width bins, 
each of which produces independent correlator output. In 
this case, there was only one bin. In this way, assuming 
the entire pulse falls within the gate, the signal-to-noise 
ratio (S/N) of the pulsar detections can be increased by 
a factor of about P/w , where P is the pulse period 
and w is the width of the gat e. The pulsar eph emerides 
were derived with TEMPO^I dHobbs et al.l 1200(1) . Before 
full correlation, we conducted iterations of gate-fitting, 
using a full-period gate with 40 bins, to confirm that the 
pulse profile was stationary over the time-range of a pulsar 
observation within an epoch and to optimize the choice of 
the gate start/stop parameters. For these three pulsars, 
the gate widths used were typically in the range of 4 — 8% 
of a pulse period, leading to gating gain factors of ~ 3 — 5. 

We performed a mostly automated data reduction 
and calibration procedure relying on the NRAO Astronom- 


1 http://www.atnf.csiro.au/research/pulsar/tempo2/ 


ical Image Processing S ystem (AI PS) a nd the scripting 
language ParselTongue dKettenis et al.l 120061) . Removal 
of data affected by radio frequency interference (RFI) 
was made runn i ng th e RFI-mitigation software SERPent 
dPeck fe Fenechl l2013i l and some further manual flagging. 
We first applied the system temperature and gain curve 
corrections as determined by the EVN pipeline^ and also 
corrected for the parallactic angle using the AIPS task 
CLCOR. Next we computed ionospheric corrections in 
TECOR with the help of total electron content maps 
as published by the Center for Orbit Determination in 
Europe^. In a first calibration run, we solved for visibility 
rates, phases, and delays in FRING for all calibrator 
sources assuming a simple point source model. Next, we 
self-calibrated on each source, improving the S/N by a fac¬ 
tor of five to ten. For each calibrator we then concatenated 
the calibrated data from all four epochs. This dataset was 
imaged to produce a global model of each calibrator source. 
The dominant CLEAN components of each source model 
were then used as the input model parameters in a second 
FRING-run. In this way, we eliminated any systematics 
caused by source structure that affected the position of the 
calibrator sources in between epochs. For each of the three 
pulsars the calibration solutions of the primary calibrator 
were of much higher quality than those of the secondary 
calibrator and, hence, were applied to the respective target 
pulsar (Tabled). 

We used the calibration solutions from the secondary 
calibrator to provide independent checks on the achieved 
astrometric accuracy as a function of angular separation 
between tar get and calibrator sour ce (‘calibrator throw’), 
as done bv iGhatteriee et all (2004). Similarly to the ear¬ 
lier data, our observations imply an almost linear increase 
in astrometric accuracy with decreasing calibrator throw 
(Fig. P- We fitted a power law, y = a ■ x b , to all data 
points obtained from 5 GHz observations, which yielded 
(a, b) = (0.15 ± 0.09, 1.17 ± 0.31). This relation, however, 
seems to only hold for angular separations of up to four 
degrees between calibrator and target. 

3. Estimates for astrometric parameters of 
B1929+10, B2020+28, and B2021+51 

We measured the position of the pulsars in each epoch 
by fitting a 2D Gaussian to the brightness distribution 
in the image plane (Fig. P using the AIPS task IMFIT. 
As a result of the high S/N (~ 150) and small beam 
size (9 ~ 3 x 5 mas), the formal errors are very small 
(9/(2* S/N) ~ 10 x 20/ras), certainly underestimating the 
real positional uncertainties. In addition to these random 
errors, residual systematic errors caused by the calibrator 
throw (Fig. P, for instance, need to be taken into account. 
A good estimate for these systematic errors is the decon¬ 
volved size 9 d of the pulsar, which is zer o for a true point 
source . Following the scheme reported in IChatteriee et HI 
(1200111 . we estimated the systematic uncertainties using 
the quantity 9 d / yj (N ant - 1) * t ohs /ti ono , where N ant = 11 
is the typical number of antennas involved, t 0 b s = 110 min 
is the total amount of time spent on each pulsar, and 
^iono = 6 min is the empirically determined atmospheric 


2 http://www.evlbi.org/pipeline/user_expts.html 

3 ftp://ftp.unibe.ch/aiub/CODE/ 
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Fig. 2: Greyscale plots of the fitted pulsar images. Rows from top to bottom show B1929+10, B2020+28, and B2021+51. 
Columns from left to right correspond to epochs one to four. Overlaid contours increase in steps of 20 percent of the 
peak flux density, where negative values are indicated by dashed contours. The absolute flux density scale (mJybeam -1 ) 
is indicated above each individual panel, and the beam size and position angle are indicated in the bottom right. The 
larger beam sizes for B2020+28 and B2021+51 compared to that of B1929+10 are due to flagging of baselines to Sh 
and/or Ur. 


Table 3: Measured positions at MJD 55629. 


Pulsar 

RA (J2000) 

Dec (J2000) 

S/N 

beam size [masXmas] 

B1929+10 

19:32:14.021289(1) 

10:59:32.90137(5) 

155 

1.01 x 5.94 

B2020+28 

20:22:37.06758(1) 

28:54:22.7563(2) 

142 

3.76 x 7.36 

B2021+51 

20:22:49.85890(1) 

51:54:50.5005(1) 

156 

1.72 x 3.97 


Notes. Numbers in brackets indicate the uncertainty in the last digit. The larger beam size for B2020+28 is due to flagging of 
baselines to Sh and Ur. 


coherence time at 5 GHz. For the total positional uncer¬ 
tainty we added both the formal and the systematic errors 
in quadrature. Table [3] lists the measured positions of all 
three pulsars in the third epoch at MJD 55629. 

To estimate each pulsar’s proper motion and paral¬ 
lax, we performed a weighted least-squares-fit to the 
measured positions. Here, we measured both parameters 
in three ways: we considered our position measurements 
alone (Fig. O, we combined our data with those of the 
publications listed in Table H] (Fig. H|), and we employed 
a bootstrapping technique. For the latter, we randomly 
sampled the position measurements that are available 
for each individual pulsar. These positions were then 


fitted and the results were stored. This procedure was 
repeated 10 5 times, yielding distributions as shown in 
Fig. [5] During the fitting procedure we allow for absolute 
positional offsets between the different data sets (typically 
of the order of several mas). Such offsets are expected 
for several reasons: i) the observations were conducted 
at different frequencies; ii) the different campaigns used 
different calibrator so urces, which may have different 
systematic errors (e.g. iKovalev et al.l 120081 : iPorca IHOOI; 
ISokolovskv et al.l 1201 lD in their ti es to the Inter national 
Celestial Reference Frame fICRF2. lMa et ahll2009l) : iii) the 
data were obtained at times that are up to ten years apart 
during which improvements to the correlator models and 
Earth orientation parameters introduce offsets; and iv) the 
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Table 4: Previous proper motion and parallax estimates and derived values. 




US 

7r 

d 

Hi 


Pulsar 

[mas yr -1 ] 

[masyr -1 | 

[mas] 

[kpc] 

[kms -1 ] 

Reference 

B1929+10 

99.0 ± 12.0 

39.0 ±8.0 

4.00 ±2.0 

n 24+ 0 09 
u<z ^-0.12 

124 + _\f 

1 

B1929+10 

94.82 ±0.26 

43.04 ±0.15 

3.02 ±0.09 

n qq+ 0.01 
u.oo_ 0> oi 

1621° 

2 

B1929+10 

94.09 ± 0.11 

42.99 ±0.16 

2.77 ±0.07 

0 361 +0 ' 009 
u.oui_ 0 o 0 9 

176±| 

3 

B2020+28 

-4.38 ±0.53 

-23.59 ±0.26 

0.37 ±0.12 

2 7Q+0.64 
z - ‘ u -0.96 

256l^ 5 

2 

B2021+51 

-5.23 ±0.17 

11.54 ±0.28 

0.50 ±0.07 

2 00+°- 23 
z<uu -0.31 

149121 

2 


References. Ill iHoogerwerf et al.l (120011 1: f 21 iBrisken et all II2002H : f 31 IChatteriee et al.l (1200411 


Table 5: Astrometric results and derived values from the estimation strategies. 



data 

sets® 


Ha 


7 r 


d 

Hi 

Pulsar 

F b 

B c 

Nobs 

[masyr -1 | 

[masyr -1 | 

[mas] 

xL 

[kpc] 

[kms -1 ] 

B1929±10 

C 


4 

94.11 ±0.52 

43.64 ± 1.06 

2.79 ±0.14 

1.34 

0 338+ 0,017 
U.OOO_o 019 

17511S 


C c , C 


10 

94.04 ±0.12 

43.39 ±0.23 

2.79 ±0.08 

1.03 

0 338+ 0 010 
U.OOO _ 0 010 

175l| 


L, C 


9 

94.20 ±0.20 

42.93 ±0.26 

2.78 ±0.05 

0.51 

0 360 + ° 006 
U.OUU_o 006 

175!^ 


L, C c , C 


15 

94.06 ±0.09 

43.24 ±0.17 

2.78 ±0.06 

0.77 

0 360 + ° 007 
U.OOU_o 008 

176!" 



C c , C 

10 

94.07lo.2o 

43 41 + 0.10 
^± 0.^1 — 0.12 

z - 1 °-0.14 


0 362 +0 019 
u.ouz _ 0 012 

177ly 



L, C 

9 

94.23lo.24 

42.95!°;|* 

2 7Q+0 - 07 
z - ‘ y - 0.12 


0 362 +0 019 
u.ouz _ 0 012 

1751® 



L, C c , C 

15 

94 .O 8 I 0 I 7 

43.251q.i3 

9 77 +O.O 8 
z - 1 1 -0.07 


0 361 +0 009 
u.oo±_ 0 .oio 

1761° 

B2020±28 

C 


4 

-3.34 ±0.05 

-23.65 ±0.11 

0.72 ±0.03 

0.05 

1 qq+0.05 

l.oy_006 

1341° 


L, C 


9 

-3.46 ±0.17 

-23.73 ±0.21 

0.61 ±0.08 

0.53 

1 63+ 0 - 18 
l.UO_0.23 

158120 



L, C 

9 

_Q 4 C+O.I 6 

-23.70lg;ll 

0.60!°+ 


1 66 +0 ' 50 
l.UO_029 

160118 

B2021±51 

C 


4 

-5.08 ±0.42 

10.84 ±0.25 

0.80 ± 0.11 

1.54 

1 9 r+0.14 
i - zo -0.17 

87l 15 2 


L, C 


9 

-5.03 ±0.27 

10.96 ±0.17 

0.78 ±0.07 

0.90 

1 OQ+ 0.10 

■*■•^- 0.12 

901 11 



L, C 

9 

— ^ 01 +°- 17 
°- ul - 0.20 

io.99i8;ii 

n 77 +O.O 8 
u -' '- 0.10 


1 30+ 019 
-L.OU_ 0.13 

9ll}° 


Notes. H) C refers to the measurements obtained in this campaign, Cc denotes the data set from Chatterjee et al. (2004), and L 
indicates that the measurements from Brisken et al. (2002) were included in the analysis. ^ Results derived from a least-squares 
fit of the measured data. H) Median values and 68% confidence interval from fitting 10 s bootstrapped realizations of the data. 


data were taken with different instruments (VLBA and 
EVN) that use different hardware/software correlators. 
Table [5] summarizes the estimates of /i and 7 r from the 
individual fits, from the different combinations of data 
sets, and from the bootstrapping method (where we quote 
the most compact 68 % confidence interval), as well as 
the implied pulsar distances and transverse velocities. 
The latter are corrected for solar motion and differential 
Galactic rotation and refer to the local standard of rest 
(LSR). Regardless of estimation strategy and combination 
of available data, all measured values are consistent within 
their uncertainties at the one-sigma level. 

For pulsar B1929+10 our re sults confirm the measure¬ 
ments of lChatteriee et all (120041 ), especially in combination 
with the earlier data. For pulsars B2020+28 and B2021+51 
our observations indicate that they are located at a distance 
of 1.39 ^)q'qq kpc and 1.25lo+kpc. Hence, they are about 
1.1 and 0.7 kpc closer to the solar system than what was 


implied by the measurements of lBrisken et all (2002) alone 
(Table [4j). Considering this discrepancy of a factor of about 
2 , we suspect that the uncertainties on the position mea¬ 
surements for B2020+28 and B2021+51 as reported by the 
authors were underestimated. Accordingly, we did not in¬ 
clude these data in the further analysis. Thus, in the follow¬ 
ing, for B1929+10 we adopted the astrometric parameters 
obtained from bootstrapping all available data, while for 
B2020+28 and B2021+51 we used the parameters as mea¬ 
sured with our new 5 GHz data alone. Hence, the analysis 
below is based on the astrometry as listed in Table [ 6 ] 

4. Simulations of pulsar orbits 

To shed new light on possible common origins of 
B1929+10/C Oph and B2020+28/B2021+51, we used the 
pulsar astrometric parameters described above and traced 
their orbits back in time through the Galactic potential. 
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Fig. 3: Relative measured positions of B1929+10 (left), B2020+28 (middle), and B2021+51 (right) with the best-fit 
proper motion removed. The solid line is the best-fit parallax from our EVN 5 GHz observations. 




Fig. 4: Same as Fig.[3j but for the fits perfor med usin g all available data. Solid dots show the most recent measurements at 
5 GHz , open stars are the data taken from iBrisken et ahl (2002), and squares are the data adopted fronr lGhatteriee et al.l 
(120041) . For better illustration, we omitted the time span without any observations. 


For the runaway star ( Oph we used the latest proper 
motio n and parallax measurements from Ivan Leeuwenl 
(120071) (Table © and adopted the va lue for the radia l 
veloci ty R a d = —9.0 ± 5.5kms _1 from Kharch enko et ai] 
( 2007 '!. Our astrometric measurements yield information 
about the transverse motion of the pulsars, but they do 


Table 6: Astrometric parameters used in the simulations 


Source 

[masyr -1 ] 

fL5 

[rnasyr^ 1 ] 

7r 

[mas] 

B1929+10 

94.08 ±0.17 

43.25 ±0.16 

2.77 ±0.08 

C Oph° 

15.26 ±0.26 

24.79 ±0.22 

8.91 ±0.20 

B1952±29 b 

-30 ±6 

-34 ±8 

1.4 ± 1.0 

B2020+28 

-3.34 ±0.05 

-23.65 ±0.11 

0.72 ±0.03 

B2021+51 

-5.08 ±0.42 

10.84 ±0.25 

0.80 ±0.11 


Notes. O) From van Leeuwen (2007) ^ Proper motion from 
Hobbs et al. (2004), parallax from the ATNF Pulsar Catalogue 
(Manchester et al., 2005). 


not contain any information about the radial velocity. 
To estimate the full 3D velocity vector, we simulated 
the possible radial component from our measured trans¬ 
verse components and the space velocity distribution 
of you ng pulsars as empirically derived by iHobbs et al.l 
(2005 ). To account for the uncertainties of the measured 
parameters p a , ps, 7r, and the unknown radial velocity, 
we assumed that all parameters are distributed normally 
(where the half-width of the Gaussian is given by the 
higher absolute value of the upper and lower errors of 
the bootstrapping results) and performed three million 
Monte Carlo simulations. The obtained velocity vectors 
were corrected for the solar motion with respect to the 
LSR, for differential Galactic rotation, and also for the 
velocity of the LSR. The Galactic potential we used in our 
si mulations is the potentia l that was described in full detail 
in lVlemmings etahl (20041) . the main parameters of which 
we summarize here in brief. For consistency reasons, the 
pulsar orbits w ere traced back th r ough the same Stackel 
potential as in lYlemmings et all (12004 ): consisting of a 
thin disk, a thick disk, and a halo component whose axis 
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Fig. 5: Probability density functions of fitting results from bootstrapping all available data. From top to bottom: 
B1929+10, B2020+28, and B2021+51. Columns are from left to right: proper motion in RA, proper motion in Dec, 
and parallax. The solid and dashed vertical lines indicate the median and the most compact 68% confidence intervals as 
listed in Table [5] 


ratios are 75.0, 1.5, and 1.02, respectively. We kept the 
relative contributions of thin and thick disks and of the 
halo at 0.13, 0.01, and 1.0, respectively. For a complete 
description of each parameter of t he St ackel potential, 
we refer to iFamaev fc Deionghel (120031). We a dopte d 


parameters for solar motion from Schonrichl (l2012f) : 


Rq = 8.27kpc and (17, V, W) = (13.84, 12.24, 6.1) kins" 1 . 
Each object’s trajectory was sampled at time intervals of 
10 3 yr using a fourth-order Runge-Kutta numerical inte¬ 
gration method. For each time step the distances between 
the two objects under consideration were computed within 
the Galactic reference frame, and we recorded only the 
simulation input parameters of trajectories that resulted in 
a minimum distance of less than 10 pc. In addition to the 


separation between the individual objects, we also com¬ 
puted their distances to the Sun (B2020+28/B2021+51) 
and to the Upper Scorpius region (B1929+10/() Oph). 
To compute the latter, we traced the trajectory of Upper 
Scorpius bac k in time using the astr ometric values as listed 
in Table 2 of lde Zeeuw et al.l (I1999H . 

Fo r consistency che c ks we used the input parameters 
of iHoogerwerf et all (|200ll ) to compute the trajectories 
of B1929+10 and C, Oph. In total, 37521 of the three 
million sampled trajectories (1.2%) cross within 10 pc 
of each other. T h is is close to the percentage found in 
IHoogerwerf et all (|200lh : 30822 out of three million, or 
1.0%. The smallest separation we found is 0.19 pc (com- 
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Fig. 6 : Probability density functions of minimal distances (left column), time of minimal approach (middle column), 
and distance to Upper Scorpius and the Sun for B1929+10/£ Oph (upper row) and B2020+28/B2021+51 (lower row), 
respectively. The stellar pair to which the figures apply is indicated in the top of each panel. 


pared to 0.35pc). Furthermore, while iHoogerwerf et all 
(l2001h reported that in 4214 (0.14%) simulations the tra¬ 
jectories of both the pulsar and the runaway star not only 
pass within 10 pc of each other, but also pass within less 
than 10 pc of Upper Scorpius, we found that 6816 (0.23%) 
of our simulations meet these conditions. The differences 
in the results are probably due to the different set-ups 
of the Galactic potentials. We di d repro duce the general 
trend found bv IHoogerwerf et al.l (120011) . however. When 
we ran the simul ations using the s a me in put parameters 
for B1929+10 as IHoogerwerf et ah (|2001 ), but used the 
latest parameters for £ Oph from van Leeuwenl (120071 ). a 
total of 82840 (2.7%) simulated orbits cross within 10 pc, 
in only 8 of which both the pulsar and the star are less 
than 10 pc away from Upper Scorpius. 

To test how much the updated solar parameters and 
the different radial velocity distribution^ influence the 
computed trajectories, we also ran the simulations for 
B2020+28/B2021+ 51 using the input p arameters for 
p a , p$, and 7 r from IVlemmings et al.l (I2004T) (Table 0. In 
our simulations 0.14% of trajectories cross within 10 pc 
(minimal distance of 0 . 10 pc), reproducing the results of 
these earlier simulations well. 


In the three million simulations that we ran using 
our bootstrapping results for B1929+10 and the latest 
astrometric parameters for £ Oph, 258272 ( 8 . 6 %) orbits 
cross within 10 pc about 0.5 Myr ago (Fig. [ 6 ]). However, 
none of these orbits yield a minimum separation of less 
than 2.4 pc, and neither the pulsar nor the runaway star 
approach the centre of Upper Scorpius to within less than 


4 We used a o ne-co mponent velocity distribution, while 
IVlemmings et al . (12004) used a two-component distribution. 


17 pc. The median radial velocity of B1929+10 required for 
it to approach £ Oph within 10 pc is 570+63 km s _1 (Fig. 
0. For completeness, we also tested the hypothesis that 
B1929+10 o nce formed a binary system with the pulsar 
B1952+29 (|Wrightl fl979l) . For th e latte r we assumed the 
proper motion from Hobbs et a l. (2004) and the parallax 
from the distance derived from the dispersion measure 
(DM) in the ATNF Pulsar Cataloguj^ and the Galactic 
electron density model from ICordes fe Laziol 120031 (Table 
0 . We assumed a parallax uncertainty of 1 mas in lieu 
of a formal error estimate in the DM-based distance. 
With these parameters, none of the simulated orbits 
crosses within 10 pc. The same is true for simulations ran 
with the same proper motion parameters, but with the 
distance estimate d = 0.42 kpc, based on the same DM 
but using instead the G alactic electron density model from 
iTavlor fc Cordesl (119931) . 


For the putative pulsar pair B2020+28/B2021+51 our new 
measurements imply a minimum possible separation of 
1.9 pc. Of the three million trajectories, 1866 (0.06%) cross 
within 10pc about I.I 6 +Q +7 Myr ago (Fig. [ 6 ]). The implied 
median radial velocities are 643+i6g f° r B2020+28 and 
433^193 km s _1 for B2021+51 (Fig.0. 


5. Discussion 

5.1. Binary companion of B1929+10 

Our new astrometric results for the pulsar B1929+10 
confirm the measurements of earlier VLBI campaigns 

5 http://www.atnf.csiro.au/research/pulsar/psrcat/ 
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Fig. 7: Probability density function (grey histograms) of astrometric parameters and required radial velocities that result 
in a minimum separation of less than lOpc between B1929+10/C Oph and B2020+28/B2021+51. Columns from left 
to right show the results for n a , 7r,and Vj a d- For the measured parameters fi a , /is, and 7r, the solid lines indicate 
the input parameter distributions derived from assuming Gaussian errors. The solid line in the last column indicates 
the input distribution for U-ad as derived from our measured transverse velocity and the empirically determined space 
velocity distribution derived bv lHobbs et al.1 (l2005h . Objects from top to bottom are B1929+10, £ Oph, B2020+28, and 
B2021+51. 


(|Brisken et al] l2002t IChatteriee et al.l l2004f) , and in com¬ 
bination with the previous position measurements, we 
place robust constraints on the uncertainties of the proper 
motion parameters. Accordingly, for our adopted parallax 
7T = 2.77jiQ Qg our values for the proper motion differ by 
more than lOcr from the parameter space that implies a 
common origin of B1929+10 and f Oph in Upper Scorpius 


in [Hoogerwerf et alJ (|200ll ). Given the new astrometry for 
the pulsar and also the updated astrometric parameters for 
the runaway star, the minimal possible separation of 2.4 pc 
between the two stars is too large to be consistent with a 
common origin of both. Moreover, the closest approach of 
roughly 17 pc to Upper Scorpius of either of the two objects 
in all of the trajectories crossing within 10 pc makes it 
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very unlikely that this region is the place of common origin. 

However, the fraction of simulated orbits that cross 
within 10 pc (~ 8.6%) is surprisingly high and implies that 
the orbits of both objects may have crossed within that 
distance about 0.5 Myr ago. The allowed range in radial 
velocities for them to pass close by, on the other hand, is 
very small and points to a very strong kick imparted to the 
pulsar at birth. Only a direct measurement of the pulsar’s 
radial velocity will further constrain the distance of closest 
approach of B1929+10 and Q Oph. 

Our data in combination with the updated astrome¬ 
try for B1952+29 make a binary origin of B1929+10 and 
B1952+29 highly implausible. 

5.2. Common origin of B2020+28 and B2021+51 


Table 7: Parallax- vs. DM-based distances 


Source 

DM 

[pc cm -3 ] 

<tr° 

[kpc] 

dNE2001 6 

[kpc] 

^TC93 C 

[kpc[ 

B1929+10 

B2020+28 

B2021+51 

3.180 (4) 
24.640 (3) 
22.648 (6) 

0 361 +0 ' 009 
u.ooi_ 0 010 

1 QQ+0.05 

1 ■°^_0.06 

1 - zo -0.17 

0.34 

2.11 

1.94 

0.17 

1.30 

1.22 


Notes. O) Parallax-based distances derived in this work. ^ Dis¬ 
tance estimate based on the Galactic electron density model 
from Cordes and Lazio (2003) O) Distance estimate based on 
the Galactic electron density model from Taylor et al. (1993). 
Except for d n , all values were taken from the ATNF pulsar 
catalogue. 

the lines of sight to B1929+10, B2020+28, and B2021+51, 
respectively. 


IVlemmin gs et a h (I2004T) u sed the astrometric parameters 


Brisken et al.l ( 20021) (Table to infer a putative 


from 

common origin of B2020+28 and B2021+51. This conclu¬ 
sions seems plausible considering that the pulsars’ 2D mo¬ 
tions lie in apparent ly opposite directions i n Galactic co¬ 
ordinates (Fig. 1 in IVlemmings et alJ120041) and also be¬ 
cause the pulsars have very similar characteristic ages of 
2.88 (B2020+28) and 2.75Myr (B2021+51). Nevertheless, 
our new proper motion measurements, in conjunction with 
our parallax measurements, which place both pulsars al¬ 
most twice as close as previous distance e stimates, rule out 
a com mon origin for these two objects. IVlemmings et all 
(|2004l) determined the percentage of orbits crossing within 
10 pc for a known binary disrupted 1 Myr ago as a function 
of astrometric uncertainties (see their Fig. 2). These models 
indicate that our improved errors should have yielded 1% 
of crossing orbit realization (within 10 pc) for B2020+28 
and B2021+51. However, in our simulations only 0.06% of 
trajectories cross within that distance, and none of the or¬ 
bits yield an approach of less than 1.9 pc. Even if we use 
the bootstrapping results with their larger errors from Ta¬ 
ble [5j only 0.08% of the orbits cross within 10 pc. Further¬ 
more, the orbits approaching each other within 10 pc do 
so at a median distance of 0.64lg’ ( j )! j ) kpc to the solar sys¬ 
tem. Given the e stimated extent of the C ygnus Superbubble 
of 0.7 — 2.5 kpc (|Vlemmings et al~]|2004l) . a common origin 
within this region is ruled out. 


5.3. Implications for Galactic electron density along the lines 
of sight 

In Table [7] we list the distances inferred from our par¬ 
allax measurements in comparison with those implied by 
the DM and the Ga lactic electron density mo d els of 
I Cordes fc Laziol (1 200. 'll . NE2001) and iTavlor et al.l (|l993l 
TC93, both obtained from the ATNF pulsar catalogue). 
While the NE2001-distance agrees with our measurement 
for B1929+10, the same model overestimates the distances 
to both B2020+28 and B2021+51 by a factor of about 1.5. 
The distance estimates for the latter two pulsars as given 
by the preceding model, TC93, agree well with our results, 
however. Hence, in combination with the DM as listed in 
Table [TJ our parallax measurements imply a mean elec¬ 
tron density of 8 . 8 (^ 2 , 17.7t(X6, and 18.lj^;|cm -3 along 


6. Conclusions 

Based on our new astrometry for pulsars B1929+10, 
B2020+28, and B2021+51 obtained with the EVN at 
5 GHz, we rule out previously proposed common origin sce¬ 
narios for all three sources. Our Monte Carlo simulations 
of the past trajectory of B1929+10 throughout the Galac¬ 
tic potential show now indication for the pulsar to have 
once been in a binary system with the runaway star () Oph 
in Upper Scorpius. Similar simulations for B2020+28 and 
B2021+51 also rule out a binary origin of the pulsars in the 
Cygnus Superbubble. 
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